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Abstract 

The mixed phase of extremely type-II layered superconductors in perpendicular magnetic 
field is studied theoretically via the layered XY model with uniform frustration. A par- 
tial duality analysis is carried out in the weak-coupling limit. It consistently accounts 
for both intra-layer (pancake) and inter-layer (Josephson) vortex excitations. The main 
conclusion reached is that dislocations of the two-dimensional (2D) vortex lattices within 
layers drive a unique second-order melting transition at high perpendicular fields between 
a low-temperature superconducting phase that displays a Josephson effect and a high- 
temperature "normal" phase that displays no Josephson effect. The former state is best 
described by weakly coupled 2D vortex lattices, while the latter state is best characterized 
by a decoupled vortex liquid. It is further argued on the basis of the duality analysis that 
the second-order melting transition converts itself into a first-order one as the perpendicu- 
lar field is lowered and approaches the dimensional cross-over scale. The resulting critical 
endpoint potentially accounts for the same phenomenon that is observed in the mixed 
phase of clean high-temperature superconductors. 

PACS Indices: 74.60.-w, 74.25. Dw, 74.25. Ha, 74.60.Ge 
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I. Introduction 



The nature of the thermodynamic phase diagram in clean high-temperature super- 
conductors in external magnetic field has been elucidated only recently. 1 Experimental 
evidence for a first-order melting transition of the Abrikosov vortex lattice has been ob- 
tained from the measurement of a jump in the magnetization and from the measurement of 
a latent heat. 2 ' 3 These two quantities are found to satisfy the Clausius-Clapeyron relation, 
which indicates that the first-order transition is indeed thermodynamic. Although the ex- 
istence of such a melting transition was expected on theoretical grounds due to the extreme 
type-II and anisotropic character of high-temperature superconductors, 4 the multi-critical 
phenomenon exhibited by the associated phase diagram for the mixed phase was not. In 
particular, the first-order melting line begins at (or near) the zero-field transition, but it 
ends strangely in the middle of the field versus temperature (T - H±) plane. The critical 
end point, furthermore, coincides (or is near) a multi-critical point, from which a nearly 
field-independent "vertical" depinning line of the vortex lattice begins at higher fields. 5 ' 6 
A weakly temperature-dependent "horizontal" second-peak line also continues on from the 
multi-critical point to lower temperatures. Bulk pinning becomes appreciably stronger at 
fields above this line. 7 

The first-order melting observed in clean high-temperature superconductors is com- 
monly interpreted theoretically in terms of elastic medium descriptions of the vortex 
lattice. 4 The melting point in such theories is usually determined by the Lindemann cri- 
terion, however, which cannot predict the nature of the transition. In other words, the 
theory is compatible with either a first-order, a second-order, or even a cross-over tran- 
sition. Elastic medium theory has also been used to describe the mixed phase of layered 
superconductors in particular. 8 A first-order decoupling transition is predicted to occur 
in the extreme type-II limit, beyond which absolutely no Josephson coupling remains 
in between the layers. 9 A subsequent high-temperature analysis of the anisotropic XY 
model with uniform frustration has demonstrated, however, that a small degree of local 
Josephson coupling must survive the decoupling transition. 10 ' 11 Last, it is important to 
remark that most elastic medium calculations to date have been performed in the gaussian 
approximation. 4 They therefore neglect important topological excitations of the vortex 
lattice, such as dislocations, that can drive the melting transition. 12 ' 13 
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Other more numerical approaches to the theoretical description of the mixed phase 
in clean high-temperature superconductors have succeeded in obtaining a clear first-order 
melting transition. These include Monte-Carlo simulations of both the frustrated XY 
model 14 ' 15 and of Ginzburg-Landau theory. 16 The former simulations are plagued, how- 
ever, by equilibration problems at low temperatures due to the artificial pinning caused by 
the underlying XY model lattice. 11 The latter simulations of Ginzburg-Landau theory, on 
the other hand, are performed in the lowest-Landau-level approximation. This approxima- 
tion is in principle valid only in the vicinity of the mean-field H c i transition, which remains 
far from the vortex-lattice melting line observed experimentally in high-temperature su- 
perconductors because of anisotropy. 2 ' 3 

In this paper, we shall develop a theory for decoupling in the mixed phase of ex- 
treme type-II layered superconductors that is based on an analysis of the corresponding 
layered XY model with uniform frustration. The analysis is effected through a partial 
duality transformation that leads to a useful neutral Coulomb-gas (CG) description for 
the Josephson coupling. 17-20 Unlike some of the Coulomb-gas analyzes of the layered XY 
model that have been performed in the past, 21-24 both gaussian and topological excita- 
tions of the vortex lattice within layers are consistently accounted for in this approach. 
This consistency is responsible for the primary result of the paper, which is that there 
exists a unique second-order transition in the weak-coupling limit that separates a cou- 
pled superconducting phase at low temperature from a decoupled "normal" state at high 
temperature. In the context of the mixed phase, the weak-coupling limit is reached at 
field components perpendicular to the layers that are large compared to the dimensional 
cross-over scale B\ = $0/^0 ( see Fig- 1)- Here, A denotes the Josephson penetration 
length. Note that the criterion used throughout to discriminate between a coupled and a 
decoupled state is the presence or absence of a macroscopic Josephson effect, respectively. 
The above result implies that neither the Friedel scenario, 4 ' 25 which is a state composed of 
decoupled superconducting layers, nor a "line-liquid" state, 14 ' 24 ' 26 which we define here as 
a set of normal layers in external magnetic field that exhibit a Josephson effect, are likely 
to be thermodynamic states in the mixed phase of extreme type-II layered superconduc- 
tors when disorder is absent. The prediction of a unique phase transition in the uniformly 
frustrated XY model with anisotropy is notably consistent with recent extensive Monte 
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Carlo simulations. 15 It was also anticipated by various authors in studies of the layered 
XY model without perpendicular frustration. 21 ' 22 ' 27 ' 28 

The duality analysis mentioned above is also employed in the paper to map out the 
phase diagram of the layered XY model with uniform frustration. A finite number of 
layers is assumed. The weak-coupling limit at high perpendicular fields is found to contain 
a low-temperature phase made up of independent two-dimensional (2D) vortex lattices 
that sustain a Josephson effect. This phase is separated from a decoupled liquid of planar 
vortices at high temperatures by a second-order melting transition akin to that shown by 
an isolated 2D vortex lattice. 29-33 The existence of such a second-order melting transition 
is demonstrated by making note of a formal equivalence between the Coulomb gas descrip- 
tions of the layered XY model with and without frustration. 20 ' 34 Also, the former smectic 
2D vortex-lattice phase is likely to be a type of "supersolid" matter (see sect. IV. C and 
ref. 13). The above weak-coupling analysis is then shown to break down at sufficiently low 
perpendicular fields of order the dimensional cross-over scale B]_. It is argued that this 
breakdown signals a cross-over transition to a flux-line lattice regime at low fields and tem- 
peratures, while that it signals a first-order decoupling transition to a potentially defective 
flux-line lattice state at low fields but high temperatures. It is further argued on this basis 
that the first-order decoupling transition terminates at a critical endpoint at a temper- 
ature of order the 2D vortex-lattice melting temperature and at a perpendicular field of 
order many times the dimensional cross-over scale B]_. The resulting phase diagram (Fig. 
1) is strikingly similar to those reported for the mixed phase of clean high-temperature 
superconductors. 2 ' 3 ' 5-7 

The paper is organized in the following way. Section II contains the weak-coupling 
duality analysis of the uniformly frustrated XY model. The principal input is the nature 
of the phase coherence in isolated layers, which is assumed to be either quasi long range 
[Eq. (17)] or short range [Eq. (18)]. The output of the weak-coupling analysis is that the 
former implies the presence of a Josephson effect, whereas the latter implies its absence. 
Equations (43) and (55) for the local Josephson coupling in the respective coupled and 
decoupled phases summarize these results. This theory is then used to determine the phase 
diagram of extremely type-II layered superconductors in perpendicular magnetic field in 
section III. The principal new result is the prediction of a coupled 2D vortex lattice phase 

4 



at low temperatures and high fields that notably displays a Josephson plasma resonance 
with exponential temperature dependence [see Fig. 1 and Eq. (66)]. A comparison of these 
results with previous experimental and theoretical work is made in section IV, while some 
conclusions are reached in section V. Finally, technical issues that concern the long-range 
nature of phase correlations in the 2D XY model and that concern a fermion analogy 
employed to describe the Josephson coupling between layers are treated in Appendix A 
and B, respectively. 

II. Theory 

The thermodynamics in the interior of the mixed phase of extremely type-II super- 
conductors (Al — > oo) can be described at least qualitatively by the uniformly frustrated 
XY model over the cubic lattice. 14 ' 15 ' 2 In the case where the superconductor is composed 
of N Josephson coupled layers, the corresponding kinetic energy for the superfluid flow 
reads 

N 

E xy = ~J\\ EEE cos[A M 0(r, /) - A»(r, /)] 

1=1 r H=x,y 
N-l 

i=i r 

Above, <f>(r, I) is the superconducting phase at a point r in layer I, while = (0, b±x, —b\\x) 
is the vector potential. Also, denotes the nearest-neighbor difference operator along 
the fi direction of the square lattice. The magnetic induction parallel and perpendicular to 
the layers is related to the frustration b through the respective identities B\\ = ($ /27r<i)&|| 
and B± = ($ /27ra)6_i_. Here the length scales d and a denote respectively the sep- 
aration between adjacent layers and the lattice constant for each square-lattice layer. 
The latter is of order the superconducting coherence length at zero temperature. Also, 
J|| = (<i/47r)($o/27rAL) 2 is the local phase rigidity within layers, while J±/a 2 denotes the 
local Josephson coupling energy per unit area. It is important to observe that the Joseph- 
son penetration length at zero temperature, Y Q a, provides a natural scale for the XY 
model (1) in the limit of weak coupling between layers, in which case the model anisotropy 
parameter Y = (J\\/J±) 1 ^ 2 is much larger than unity. 4 This shall be assumed throughout. 
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Last, we remind the reader that any generalized phase auto-correlation function set by an 
integer source field, p(r), is related to the corresponding partition function 

Z f Y \p] = J 2ty e -*&/*BV5> (2) 

by the quotient 17 ' 18 

; e xp[z^p(r)^(r)]) = zf Y \p]/zf Y [% (3) 



A knowledge of these amplitudes characterizes the thermodynamics of the layered XY 
model (1). 

A. Layered Coulomb Gas. We shall now employ the well-known dual representation 
of the XY model (1) that is based on the Fourier series expansion 18 



in9 



e p cos9= I H ((3)e 



of the Gibbs distribution in terms of the modified Bessel functions I n (x). This identity 
allows the phase variables to be integrated out of the partition function at the price of 
introducing a new integer field n that lives on each link of the XY model (1). Substitution 
into Eq. (2) thereby results in the dual form 

z { * Y \p} = if(p 1 _)i* Af (p ll ) n^J^A^i.-KOj-n^^^^^e-^w^w (4) 

W(r)} v 

expressed in terms of modified Bessel functions and the quotients t n (/3) = I\ n \(f3) / I (f3). 
Above, n M (r) is an integer link-field on the layered lattice structure of points r = (r, /), 
with [i = x,y,z, and (3^ = J^/ksT. Also, Af denotes the total number of sites, while 
M' = J\f(l — iV _1 ) gives the total number of rungs between layers. [The perpendicular 
link fields at the boundary layers are set to n 2 (f, 0) = = n z (r, N)]. We now take the 
crucial step in the theory by observing that only the configurations with n z (r) = 0, ±1 
are relevant in the weak-coupling limit, 35 ' 36 f3± 1. After re summation over the parallel 
link fields n x and n y in Eq. (4), we then obtain the form 

Z ( i ) Y \p]=lfm E yo M -^Z ( p Y [ qi ]-e-^^ (5) 
{Mr)} 



for the partition function of the 3D XY model in terms of sums of products over partition 
functions 

Z^ Y [ qi ] = J Vfr e -*&/*B*y X>*« (6) 
for frustrated 2D XY layers in isolation, where 

E xy = ~J\\ E E cos[A^,(f) - A,{r)\ (7) 

r fJ-=x,y 

is the intra-layer superfluid kinetic energy, and where 

qi(r) = p(f, I) + n z (f, I - 1) - n z (f, I) (8) 

is the source integer field. Above also, N[n z ] = J2ri \ n z(r, /)| is the number of n z (fluxon) 
charges per configuration and yo = Ii(P±)/Io(P±) is the bare fugacity. The latter tends to 

Vo = P±/2 (9) 

in the weak-coupling limit <C 1. By analogy with Eq. (3), we now identify the quotient 
^xyfel/^xyM with the generalized auto-correlation function 

d[ qi ] = (exp[t5>(f)*(f;0]) (10) 

r 

of the 2D XY model with frustration. 17 ' 18 Comparison with Eqs. (5) and (6) thus yields 
the form 

Zfy [P\ = if (#]_) " Z CG [p] ■ nl.zfy [0] (11) 

for the partition function of the layered XY model in terms of a product of a layered 
Coulomb gas ensemble 

Z CG [P\= E yo [nz] -^liCih]-e-^ n ^ (12) 

{n z (r)} 

with iV XY model layers in isolation (see ref. 27). This is the final result of taking the 
selective high-temperature limit, yo — > 0, which is assumed throughout. 

Before proceeding further, we now derive a useful relation between the density of n z 
charges and the local Josephson coupling. Eq. (1) indicates that the latter is given by 

(cos0 w ) =Af'- 1 dlnZ^ Y [0}/dp ± , (13) 
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where 0/ 5 / + i(f) = <fi{r, l + l) — 0(r, /) — A z (r) is the gauge-invariant phase difference between 
consecutive layers. Yet the Coulomb-gas ensemble (12) yields the identity (N[n z ]) = 
yo(d\nZcG[Q]/dyo) for the total number of fluxons on average. By the factorization (11) 
for the layered XY model, we thereby obtain the general expression 

(N[n z ])/Af' = 2 Z/o «cos0 M+1 ) - y ) (14) 

for this quantity per rung in terms of the local Josephson coupling. 

To demonstrate that (12) is indeed a layered Coulomb gas ensemble, consider now 
a single neutral pair of unit n z charges that lie in between layers /' and I' + 1 in the 
absence of an external source, p = 0, with the negative and positive charges located 
at planar sites r\ and f 2 , respectively. This represents the first non-trivial term in the 
selective high-temperature expansion (12). The gauge-invariant product over intra-layer 
autocorrelation functions in the layered Coulomb gas ensemble then reduces to the product 
Ci> (fi, r 2 )Cp +1 (ri, r 2 ) of the corresponding two-point functions, 



Ci(r u f 2 ) = (exp i(p(f l7 l) -i(p(f 2 J) ) , (15) 

\ L J / J j_ — 

within isolated layers. This function takes the form 

CH^fs) = IQin-r^e-'X 1 '^ (16) 

in terms of a suitably gauge-transformed vector potential A', and in terms of a magnitude 
that varies with the separation as 

\Ci(^\=go(ro/\f\) V2D for |f| « £vx, (17) 

and as 

\C l (r)\=g eM-\r\/U) for |f| » £vx- (18) 

Here, 

V2D=Vsw+Vv X (19) 

is the 2D correlation exponent inside layer /, where ?y sw = (27r/3|| ) 1 and t] vx are the 
respective spin- wave and vortex contributions (see Appendix A). Also, £ vx denotes the 2D 
phase correlation length, while the length r = a vx /2 3 / 2 e 7E is set by the inter- vortex scale, 

a vx = ($o/5±) 1/2 , (20) 



and by Euler's constant, 18 7e- General scaling considerations 37 indicate that the ratio 
£vx/tt V x depends only on temperature. Last, the prefactor in expression (17) is of order 
go ~ (a/dyx)^- At fields B± > $ /7o 2 a 2 , it is therefore of order unity for anisotropy 
parameters, Y < e 27r ^" , that are not astronomically large. The effective layered CG 
ensemble (12) therefore takes form 19 



"'EE n,(f,/)A 2 (f,/)l (21) 

* r J 

in the limit of dilute fluxon (n z ) charges. Above, the Coulomb gas ensemble (12) has been 
coarse grained up to the natural ultra-violet scale a vx . In particular, the sums above are 
restricted to a square sublattice with lattice constant a vx . This requires the introduction 
of an effective coarse-grained fugacity 

V = 9o(a vx /a) 2 y . (22) 

At relatively small separations, \r\ — r*2| <C £ vx , within a 2D correlation length, the fluxons 
experience a pure {vj^} = 0) Coulomb interaction set by i]2D [see Eq. (17)]. At large 
separations \f\ — r*2| £ vx , on the other hand, the fluxons experience a pure (t]2D = 0) 
confining interaction (ri, f^) = \r\ — r 2 1 / ^vx [see Eq. (18) and refs. 35, 38 and 39]. 

Study of the asymptotic behavior exhibited by general n-point auto-correlation functions 
(10) for the 2D XY model reveals that the form (21) for the Coulomb gas ensemble (12) 
remains valid in the case of a general fluxon charge configuration, n 2 (r, /) (see Appendix 
A) . However, the preceding confining interaction exists only between those well separated 
points r*i and r*2 in a given layer / that are connected by a string. 

The coarse- graining (21) of the original Coulomb gas ensemble (12) assumes implic- 
itly that there be no more than one n z = ±1 charge per unit sublattice area a 2 x . This 
is equivalent to the condition {N[n z ]) /M' < a 2 /a 2 x for the density of such fluxons. Com- 
parison with expression (14) for this density then implies that the preceding condition is 
satisfied at temperatures and fields such that f3±a 2 rx /a 2 < 1. The validity of the above 
coarse-grained Coulomb-gas ensemble (21) is hence assured at high perpendicular fields 
B± > /3||<E>o/7o 2 a 2 . Below, we determine the thermodynamic nature of the two phases 
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that correspond to quasi long-range (17) and to short-range (18) phase correlations within 
isolated layers. 

B. Coupled Phase. Consider first the case (17) where quasi long-range intra-layer phase 
correlations are present: £ vx = oo and V^ ing = 0. Inter-layer n z charges in the Coulomb 
gas ensemble are screened at low temperatures, r]2D < 2. Global charge conservation is 
no longer enforced in this instance. Let us also assume the weak-coupling regime defined 
by the inequality y <C 1 to be satisfied by the effective fugacity (22) of the coarse-grained 
CG ensemble (21). Following the standard prescription, 35 ' 36 we then sum independently 
over charge configurations at each site, with the restriction to values n z = 0, ±1. An 
appropriate Hubbard- Stratonovich transformation 24 of the CG (21) in the absence of a 
source (p = 0) followed by a Debye-Hiickel type of approximation 35 ' 36 reveals that it is 
equal to the corresponding one 

Z LD [0] = / V9e- E ^/ kBT (23) 



for a renormalized Lawrence-Doniach (LD) model 4 up to a factor that depends only on 
The LD energy functional 

n i 1 N-l i 

£ld = J||E E E ^0 l ) 2 -2yk B TY, E ™s(9 l+1 - 9 t - A z ) 

1=1 H=x,y r 1 = 1 f 

(24) 

is summed over the coarse-grained sublattice. Here, 

J,, = k B T/2n V2D (25) 

is the macroscopic phase rigidity of an isolated layer. 40 (The validity of the Debye-Hiickel 
type approximation will be established a posteriori.) Taking the continuum limit of the 
lattice energy (24) yields the more familiar expression 

» r n 1 N-i 

E LD =J|| / d\ £-(W0 2 - A " 2 E cos (^+! ~ e i- A ^ ( 26 ) 
^ U=l 1=1 

for the energy functional of the LD model, with a renormalized Josephson penetration 
length 

A = a( / 3||A 7o yo) 1/2 , (27) 
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expressed in terms of the parameter f3\\ = J\\/k B T. Substitution of the bare fugacity (9) 
that corresponds to the conventional XY model yields a more familiar expression 



Ao=a(J||/s </±) 1/2 (28) 

for the renormalized Josephson scale. The above continuum description (26) is understood 
to have an ultraviolet cut off «o ~ a vx on the order of the coarse-grained sub-lattice 
constant. It is also gaussian within layers, which means that parallel Josephson vortices 
are the only topological excitations that it contains explicitly. 

The above renormalized LD model is known to be macroscopically Josephson coupled 
at temperatures below /c^T^o = 4n Jy. 21-24 ' 41 ' 42 To illustrate this fact, consider first the 
minimal double-layer case, N = 2. The zero-temperature line tension of a single Josephson 
vortex is equal to the linear energy density 

£||(0) = 2 5 / 2 J||/A (29) 

of a sine-Gordon soliton. 24 The condensation energy, on the other hand, is in general equal 
to — Gcond = ksT In Zld- At zero temperature and zero parallel field, this energy per 
vertical rung is thus given by 

-^(0)/^' = J> 2 /A 2 . (30) 

At the unbinding temperature /cbTle = 2ttJ\\ for Josephson vortex/ant i- vortex pairs, 41 ' 42 
on the other hand, the double-layer LD model can be analyzed exactly through a mapping 
to ideal spinless fermions (see Appendix B). The line-tension so obtained is given by the 
expression 

£\\(Tle) = 7rJ||a /Ao, (31) 

where cto is the natural ultraviolet length scale of order the underlying lattice constant a vx 
for LD model (26). This mapping can similarly be employed to obtain the expression 

-G { 1ATle)/^' = (J||/8)(aa /A 2 ) 2 (32) 

for the condensation energy at zero parallel field. Notice that both the line-tension and the 
condensation energy above are renormalized down with respect to their zero-temperature 
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values by factors of a /A . Finally, both the fermion analogy 42 (see Appendix B) and 
the original CG ensemble (21) indicate that the renormalized LD model (26) decouples 
at /cbT* = 47rJ||. 21-24 The fermion analogy in particular yields that both the line tension 
£\\{T) and the condensation energy G^J nd (T) vanish in an exponentially activated fashion 
as T approaches T* from below (see Appendix B). A "semi-classical" calculation 43 of 
the renormalization to the zero-temperature line tension ey that results from thermal 
wandering of the Josephson vortex is consistent with the exact results listed above at the 
three temperatures T = 0, T L e, and T*. 24 This calculation obtains a line-tension of the 
form 

£||=4J||/Aj, (33) 

where Aj is a renormalized Josephson penetration length that varies with temperature 
following 



, x [l-(Wa)]- 1 

Aj = / A ' 

r \ 2V 2 r 



(34) 



The exponent above takes the values 1, 2 and oo at the respective temperatures T = 
0,Tle, and T*. This agrees with Eqs. (29), and (31) if the ultraviolet length scales are 
identified by ro = (7r/8)ao- It is also consistent with a decoupling transition at T*. 

Given the success of the "semi-classical" result [Eqs. (33) and (34)] for the line-tension 
in the double-layer case, we propose that the condensation energy per vertical rung of the 
sine-Gordon model [Eq. (26) for A = 2] generally takes the form 

-^° ) nd /^ / = / i oJ||a 2 /2A 2 (35) 

in zero parallel field, where ho is weakly temperature (t}2d) dependent and of order unity. 
Comparison with Eqs. (30) and (32), for example, yields the assignments h (0) = 1 and 
/i (7le) = 4/-7T 2 for this prefactor. Since the Gibbs free energy of a double-layer is generally 
given by the formula 

Gcond = Gcolid + (^Lyd)^/^ (36) 
in the limit of vanishingly small parallel field, B\\, we obtain the final result 

1 



-G con( l/Af = J\\ 



-ho(a/X 3 y - 4(a/Aj)(ad|B|||/$ ) 



(37) 
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for the condensation energy in such case. We shall now make use of this result in order to 
compute the local Josephson coupling (13). Since we have the identity Zcg[0] = Zld[0] 
up to a factor independent of the factorization (11) implies that 

(cos0 IiI+ i) =y + d(-G cond /M'k B T)/d/3 ± . (38) 

Substitution of the double-layer result (37) for the Gibbs free energy above yields the final 
expression 

u-l 



(cos0 M+ i) =y + h g v\ 



2rJ 

,2 




AH / \ A / B'\ 

for the local Josephson coupling in the limit of small parallel fields, where 

^=[1-(W2)]" 1 (40) 

is the temperature dependent exponent that appears in Eq. (34), where B^ = $o/A <i 
is the crossover field above which Josephson vortices overlap, and where hi = 2 3 / 2 . The 
second term on the right-hand side of Eq. (39) agrees with the result obtained by Glazman 
and Koshelev 44 for the local Josephson coupling in layered superconductors at zero field 
(r ~ a). 

Consider next the renormalized LD model (26) in the case of an infinite number of 
layers, N — > oo . The zero-temperature condensation energy at B\\ = is clearly equal 
to that of the previous double-layer case (30). In addition, mean-field treatments of the 
fermion analogy for the renormalized LD model indicate that the Gibbs free energy G con d 
at low parallel fields is in general given by the original energy functional (26), 42 but with a 
renormalized Josephson penetration length Aj(T) > Aq. In the continuum limit, we have 
for example 



Gcond ^ G^(0) + 4 J dxdydz 



(—\ 2 + - 2 (—-^ 

dx ) \dy I ^ \dz d 



(41) 



with the effective mass anisotropy parameter 7 = Aj/gL The line-tension of a single parallel 
Josephson vortex in the present case — > 00 is then given my the known result 4 



£,l =( 7 rJ||/A J )ln(A L /rf). 
13 



(42) 



Here, the resulting logarithmically divergent integral (41) has been cut off by the London 
penetration length, Al, which is taken to be the natural infrared scale. Also, the conden- 
sation energy in the absence of parallel field is then given by the zero-temperature result, 
Eq. (30), but with the bare Josephson scale Ao replaced by Aj(T). If we now assume that 
Aj(T) has the same scaling form 44 as the analogous scale in the double-layer case, Eq. 
(34), then a repetition of the previous steps yields an expression for the local Josephson 
coupling in the presence of a small parallel field of the form 

<-W^+/o(^) , -/ 1 (g)'' / (43) 

Here r = a vx /2 3 / 2 e 7E , Aj is of order A , and the prefactor fi is of order g ln(X L /d). 
This result for the local Josephson coupling is corroborated by the observation that the 
second-term above coincides with the intra-layer auto-correlation function (17) evaluated 
at a separation r = Aj of order A = 70 • d ! 44 This suggests that the prefactor and 
the exponent for the second term have the limiting values 77 — > r/ 2 D and fo — > go at low 
temperatures, 772D <C 1, which is corroborated by the double-layer result (39). The latter 
corresponds to the assignments r\ = (77^ — -|) _1 for the effective exponent, 44 Aj = Aq/2 1 / 2 
for the effective anisotropy scale, and fi = higo/(l — \t]2d) for the prefactors (i = 0, 1). 
Notice that /o(0) = go and /o(Tle) = (8/n 2 )go in this case, which indicates that fo/go is 
close to unity at low temperatures r]2D <C 1. 

We now demonstrate the validity of the Debye-Hiickel type of approximation referred 
to earlier. In passing from the layered CG ensemble (21) to the renormalized LD model 
(24), the amplitudes l+2y cos(9i + i~6i—A z ) that appear in the partition function as a result 
of the Hubbard-Stratonovich transformation 24 are replaced by exp[2y cos(6>z + i — Q\ — A z )\ 
(see refs. 35 and 36). This is valid as long as Debye screening of the fluxon charges occurs, 
which requires many such charges within a screening cloud: 7rA 2 • {N[n z ]) jM' a 2 > 1. 
Study of the previous formulas in the double-layer case, Eqs. (14), (34), and (39), yields 
the result ^(fo/ go)P\\ for the former quantity. The Debye-Hiickel type approximation is 
therefore valid at low temperatures ksT < Jy, which as we shall see in the next section 
contains the present coupled phase. 

Last, we shall compute the perpendicular phase rigidity, pjr = -<^n(G cori d/ N) |o in the 
coupled phase. Here, A' z represents the longitudinal component of the vector potential that 
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is assumed to be equal across all of the perpendicular links. Periodic boundary conditions 
in the direction perpendicular to the layers are implicit. First, the duality transformation 
(4) yields the general expression 

ri- = ^" 1 ([E n -(* ? '0] 2 )*Br ( 44 ) 

r,l 

for the perpendicular stiffness in terms of a fluxon average over the dual ensemble (21). 
Yet the n z charges are screened in the coupled phase. Such short-range correlations among 
the fluxons suggests the approximation A/" -1 (Er z n z(r, I)] 2 } = {n 2 z ) for the average that 
appears above. But since the fluxon charges are effectively restricted to take values n z = 
0, ±1 in the selective high-temperature limit, yo — > 0, we also have approximately that 
(n 2 z ) (N[n z \)/N'. Use of Eq. (14) thus yields the (mean-field) estimate 



Ps/J± = (cos0 M+ i) -y (45) 

for the perpendicular phase rigidity of the uniformly frustrated layered XY model in 
the selective high-temperature limit of the coupled phase. Comparison with the previous 
results for the local Josephson coupling indicates that the perpendicular phase rigidity 
decreases with increasing temperature and field following the sum of the last two terms on 
the right-hand side of Eq. (43). 

In conclusion, a macroscopic Josephson effect exists at low temperatures, r]2D < 2, in 
extremely type-II layered superconductors if isolated layers display quasi-long-range order 
(17). Indeed, the analysis just completed demonstrates that this is the case even in the 
weak-coupling regime (cos 4>i,i+i) «C 1 of the coupled phase reached at high perpendicular 
fields B ± > $ /7o 2 « 2 - 

C. Decoupled Phase. Consider next the case (18) in which intra-layer correlations 
are short range: £ vx < oo. Although the coarse-grained form (21) of the Coulomb gas 
ensemble remains valid in this instance, a direct analysis starting from the original form 
for the ensemble defined by Eqs. (8) and (12) is more expeditious. Since the phase 
auto-correlation functions (3) of interest are those that probe inter-layer couplings, the 
associated source field shall be assumed to take the form 

p(f,l)=n^(f,l-l)-n^(f,l), (46) 
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where n\ (r) is a fixed fluxon "charge impurity" field. Notice that the intra-layer sources 
(8) can then be simply re-expressed in terms of the net fluxon charge distribution 

m z = n^+n z (47) 

as 

qi(r) = m z (r,l - 1) - m z (r,l). (48) 

This image of fluxon charge impurities is useful in the calculation of inter-layer autocorre- 
lators. In particular, Eqs. (3), (11) and (12) yield the identity 

(exp(t5>(r)*(r) - n z °\r)A z (r)])) = Z' CG \p}/Z CG [0} (49) 

r 

between the corresponding gauge-invariant quantity and the quotient of the modified par- 
tition function 

Z' CG [P]= £ ^•nf^afel'e 4 ^^ (50) 

{n z (r)} 

with the unmodified one (12) in the absence of a source. We shall now apply this to the 
calculation of the local Josephson coupling (e l ^ ! ,+1 ), which corresponds to fixing a unit 
fluxon "charge impurity" at some point in between layers / and I + 1. In the selective 
high-temperature limit, yo — > 0, the configurations that contribute to the numerator (50) 
of the quotient are limited to those with a net fluxon charge — N[n z ] = —1 in between 
layers / and / + 1 only. Likewise, the denominator Zcg[0] can be approximated by unity 
in such case. We thereby obtain Koshelev's formula 10 ' 11 

( e i*u+i) <x yo J d 2 rd{Q,r)Ct +1 {Q,r)e-^ x /a 2 (51) 

for the local Josephson coupling in the decoupled phase. Notice that short-range auto- 
correlations (18) yield a finite integral above, while quasi-long-range autocorrelations (17) 
yield a divergent integral for 2D correlation exponents that satisfy rj2D < 1- Last, Eq. (51) 
also implies that the first derivative d(e l ^ l ' l+1 ) /dBn of the local Josephson coupling is null 
in zero parallel field. This result survives to all orders of the selective high-temperature 
expansion in the decoupled phase. Eq. (38) in turn implies that the condensation energy 
satisfies d 2 G con d/df3±dB\\ = at B\\ = 0. We then have a null line tension for Josephson 
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vortices in the decoupled phase, since the latter yields £y = = £|||j,_=o- I n conclu- 
sion, the Josephson effect is absent in the weak-coupling limit of extremely type-II layered 
superconductors when the intra-layer phase correlations of isolated layers are short range. 

We shall now demonstrate that phase auto-correlations in between layers are short- 
range in the decoupled phase. Since the introduction of parallel magnetic field generally 
reduces inter-layer phase correlation [see Eq. (55) below], it is sufficient to consider the 
zero-field case, B\\ = 0. By the previous discussion, the gauge-invariant phase autocor- 
relator (e^' i+ ™) in between n + 1 adjacent layers is equivalent to fixing a column made 
up of n > 1 unit fluxon "charge impurities" in between layers / and / + n. The lowest 
order configurations that contribute to the numerator (50) of the CG ensemble are those 
with net fluxon charge —N[n z ] = —n distributed evenly in between layers I and l + n only. 
Equations (16), (49) and (50) therefore yield the expression 



(e 



l<Pl,l + n\ — 



= (Wa 2 rnlt n 



/ 



d 2 rv\Ci'{ri> 



+ f l+1 + ... + f l+n) (52) 



for this quantity to lowest order in the fugacity. Substitution of the Fourier representation 
5^ 2 \R) = (2n)~ 2 J d 2 qe l ^ R for the 2D 5-function above in addition to replacing the phase 
factor by unity then yields the inequality 



l<Pl,l + r 



< 



Co 



d 2 q 




n+l-i 



(y C /a 2 r, 



(53) 



where 



C„= I d 2 r\C v (r v )\e^ 



(54) 



is the Fourier transform of the intra-layer auto correlator. The latter are assumed to be 
identical for each layer. [Equation (53) is an equality in the absence of field, B± = 0, or 
for n = 1.] The prefactor in brackets above (53) typically decays polynomially with n [see 
Eq. (56) below]. We then conclude that the inter-layer auto correlator (e l ^ M+ ™) decays 
at least exponentially with n in the weak-coupling limit yo — > 0. This implies that the 
macroscopic phase rigidity, pj- , in the direction perpendicular to the layers is null in the 
decoupled phase. Notice that such a null result is consistent with expression (44) and with 
the fluxon charge neutrality that is characteristic of the decoupled phase. 
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To get a more concrete idea of the results just obtained in the selective high-temperature| 
limit, yo — > 0, for the decoupled phase, we shall now assume the form (16) and (18) for 
the short-range intra-layer auto correlator. Substitution into (51) yields the expression 

2 n -3/2 , \ 2 



(cos 01,1+1 ) = 2 



1 + 



||?vx 




gov (55) 



for the local Josephson coupling in parallel field. Notice the anti-cyclotronic 1/B± depen- 
dence and 1/T dependence in the above expression for the local Josephson coupling that 
is generally expected from Koshelev's formula (51). Notice also the quadratic dependence 
with parallel field that is consistent with a null parallel line tension, £y =0. It is inter- 
esting to remark that outside of the 2D critical regime, where £ vx ~ a vx , the result (55) 
coincides with the form (43) for the local Josephson coupling in the coupled phase. The 
corresponding exponent r\ = 2 occurs precisely at the depinning temperature, r]2D = 1 5 
for Josephson vortices in the double-layer case (39) ! 42 Last, the inter-layer autocorrelator 
satisfies Eq. (53), with Fourier components C q = 2ngo^ 2 x / (q 2 ^ 2 x + l) 3 ^ 2 for the short-range 
correlations (18). We thereby obtain the inequality 

(e l ^)<^^(2nye v Ja 2 v T (56) 

for the autocorrelator in between n + 1 layers in the weak-coupling limit of the decoupled 
phase. 



III. Mixed Phase 

We shall now employ the theory developed in the preceding section for a finite number 
N of Josephson-coupled XY layers (1) with uniform frustration 14 ' 15 ' 24 to determine the 
thermodynamic nature of the mixed phase in the corresponding extremely type-II layered 
superconductor. The Josephson energy will be assumed to follow the temperature depen- 
dence Jj_ oc T c o — T in the vicinity of the mean- field critical temperature, 45 T c o, while 
the anisotropy parameter 7q = (J\\/J±) 1 ^ 2 will be assumed to be a constant that is large 
compared to unity. 

A. Phase Boundaries. Consider first the limit of weak local Josephson coupling, 
(cos 0/,/+i) — > 0. Comparison with the results (43) and (55) indicates that this limit 
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coincides with the regime of high temperatures and high fields, T 3> Tj = J±/ks and 
B± 3> B]_ = $> /Aq. It is well known that an isolated lattice of 2D vortices (6) melts at a 
temperature 

k B TS D) = J\\/20, (57) 

above which quasi long-range positional correlation is lost. 29-31 The transition is driven 
by the unbinding of dislocation pairs and it is second-order. 32 ' 33 This is reflected by the 
analysis of the Villain model (Appendix A), where the elementary excitations are explic- 
itly interstitial/vacancy pairs with respect to the triangular vortex lattice that exists at 
zero temperature. 12 ' 13 Suppose now that the inter-layer coupling is weak enough so that 
the inequality Tj <C Tm°^ is satisfied. (This is consistent with a minimum anisotropy 
parameter 7q between four and five.) Let us also make the plausible assumption that the 
phase correlations within an isolated layer inherit the quasi long-range behavior shown by 
the vortex positions at low temperature. 29 ' 30 By the analysis of the previous section, we 
then conclude that the layers exhibit a macroscopic Josephson effect at low temperature 
T < Tm following the renormalized LD model (26), while they decouple at high tem- 
perature T > Tm D \ In the limit of high perpendicular fields, which coincides with the 
present weak-coupling limit Y — > oo at fixed field, we therefore have a low-temperature 
phase composed of N 2D vortex lattices that exhibit a Josephson effect, and that is sep- 
arated from a decoupled vortex liquid phase at high temperatures by a second-order line 
at T = Tm D ^ ■ Hence, neither the Friedel scenario 25 (decoupled superconducting layers) 
nor the "line-liquid" state 14 ' 24 ' 26 (coupled normal layers) are thermodynamically possible 
in the weak-coupling limit. 

Consider next the weak-coupling regime at small yet non-vanishing local Joseph- 
son coupling, (cos0^+i) <C 1. It is important to observe first that the selective high- 
temperature expansion for the local Josephson coupling in the decoupled phase breaks 
down at perpendicular fields of order the scale 

Bl =gl(3\\{U/a^) 2 Bl (58) 

in the absence of parallel field. The corresponding approximation (55) is of order unity in 
such case, which coincides roughly with the identification of length scales £ vx ~ A . The 
selective high-temperature expansion (43) for the coupled phase, on the other hand, breaks 
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down at much lower fields of order B\ for temperatures T ~ Tm ■ For a given temperature 
inside the decoupled phase, the confining nature of the CG (12) is most prominent at 
perpendicular fields that are much larger than . In particular, the string interaction (18) 
binds together dilute fluxon-antifluxon pairs into stable dipoles of dimension £ vx that do not 
overlap in this limit. The regime is best described physically by a decoupled vortex liquid 
with short-range correlations on the scale of £ V x- 10 At weak-coupling, (cos <j>i,i+i) <C 1, Eq. 
(55) indicates that the effective fugacity y of the CG ensemble (21) is small compared to 
unity. Eq. (55) also implies that the cross-over temperature T x for a fixed perpendicular 
field B± (58) lies inside of the 2D critical regime (£ vx 3> a vx ) in such case. A similar dilute 
CG description (21) exists for the layered XY model without frustration, in which case 
coarse-graining is absent, a vx — > a, and the effective fugacity y is replaced by the bare one, 
yo- 20 By analogy, we therefore conclude that a second-order melting transition takes place 
at a temperature T m that lies inside of the dimensional crossover window Tm < T < T x 
for a fixed perpendicular field B± 3> B\_ (see ref. 34). It is worth pointing out that the 
bound fluxon pairs begin to overlap inside of this regime [see Eq. (62) below]. 
Before we continue, it is useful first to define a decoupling contour 

(cos0 M+ i) = (cos0 M+ i)r> (59) 

in the T-B± plane for fixed parallel field, B\\ . Here, (cos (f>i,i+i) d is a constant less than but 
of order unity. 11 ' 46 Its value will be determined a posteriori below. In particular, the result 
(43) for the local Josephson coupling in the coupled phase yields (i) a nearly "vertical" 
contour line at temperatures of order the Josephson energy, ksTj = J±, in the limit of 
extremely high perpendicular fields B± compared to Bj_, and (ii) a contour line at higher 
temperatures Tj T < T 7 ^ D ' ) with a perpendicular field that increases exponentially with 
T* /T as 

tf D ~Bl(/o/(cos<^ +1 ) D ) 2 /' 7 . (60) 

In the decoupled phase at high temperatures outside of the 2D critical region (£ vx ~ a vx ), 
on the other hand, the result (55) for the local Josephson coupling indicates that the 
decoupling contour lies at a perpendicular field 



H D ~ gl(3\\By (cos (j) hi+1 ) D 
20 



(61) 



that decreases monotonically with temperature as 1/T. A similar result for the decoupling 
transition is obtained using the elastic medium description. 8 ' 9 Since this high-temperature 
contour should connect smoothly with the low-temperature one (60) in the coupled phase, 
the contour line (59) must cross the 2D melting line T = Tm at a perpendicular field 
many times B±. 

Suppose now that the local Josephson coupling, (cos^^+i), begins to approach unity, 
which can be achieved by either lowering the perpendicular field or the temperature. The 
selective high-temperature expansion in the coupled phase (43) breaks down at tempera- 
tures and fields in the vicinity of the decoupling contour (59). Yet the CG ensemble (21) is 
screened in the Josephson-coupled phase at low temperatures T < for small effective 

fugacity y. The neutral plasma of fluxons is dilute in this instance by Eq. (14). Since an 
increase in the fluxon density with respect to the dilute limit can only increase the effect of 
screening, then no thermodynamic phase transition is possible as a function of field (y) for 
T < Tm D \ This general argument is corroborated by the phase diagram for the neutral 
2D Coulomb gas. 47 The break-down of the selective high-temperature expansion in the 
coupled phase that occurs in the vicinity of the decoupling contour (59) must therefore 
signal a crossover into a flux-line lattice regime. 

On the other hand, what happens in the decoupled phase as the perpendicular field is 
lowered through the cross-over field at temperatures and fields that are beyond the 2D 
critical regime (£ vx ~ a vx )? It is useful in this case to determine first the point at which 
bound fluxon-antifluxon pairs begin to overlap. Since the average fluxon occupation per 
vertical rung is equal to 2a 2 /r^, where r s denotes the average spacing in between dipoles, 
and since the second term in expression (14) for this quantity is negligible when the local 
Josephson coupling is of order unity, we then have 

tlJrl = 2/o(evx/a 2 )(cos0 M+ i). (62) 

Notice that the prefactor on the right-hand side above is of order the high-temperature 
approximation (55) for the local Josephson coupling. Since both the approximate (55) 
and the true value for (cos 4>i,i+i) are °f order unity 11 ' 46 along the decoupling contour 
(59) at temperatures outside of the 2D critical regime, Eq. (62) then indicates that the 
fluxon pairs begin to overlap (r s ~ £ vx ) in such case. The dipoles disassociate, however, 
once they overlap due to the ineffectiveness of the string (18). The system must therefore 
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eventually experience a (inverted) phase transition into a screened CG above a critical 
coupling (cos^z-i-i)^) less than but of order unity. 11 ' 46 Also, the phase transition must be 
first-order since there is no diverging length scale nearby (£ vx ~ a vx ). 48 This is consistent 
with elastic medium descriptions of the vortex lattice in layered superconductors, 9 which 
also predict a first-order decoupling transition at a perpendicular field Ho of order (61). 
Last, continuity with the second-order melting/decoupling transition that takes place at 
higher fields implies that this first-order decoupling line must end where the former line 
begins. We therefore predict a critical endpoint at a temperature and field of order Tm D ^ 
and many times B]_, respectively. 

Finally, the layered XY model (1) shows a second-order phase transition in the absence 
of frustration at a relatively large critical temperature koT c ~ Jy in comparison to the 2D 
melting temperature (57). 34 The above first-order decoupling line must therefore end in the 
vicinity of this zero-field critical point as temperature increases. Also, the previous duality 
analysis can be repeated in its entirety for the zero-field case, where the replacement a vx — > 
a must be made globally. 20 Eqs. (43) and (55) then imply that the local Josephson coupling 
is of order unity at temperatures between Tj = J±/ko and min [27rJ||/A; B ln7^, T^ 2D \ We 
therefore have To < T c , which is also consistent with the second critical endpoint in the 
vicinity of T c at zero field. 

The above results are summarized by the phase diagram shown in Fig. 1. The 
phenomenology J± = Ejq(T c q — T)/T c q for the Josephson energy in the vicinity of 
the zero-field transition at T c q yields the linear temperature dependence 2 ' 8 ' 9 Ho(T) = 
~f2 2 H c2 (T) for the first-order decoupling field (61), where H c2 (T) ~ (<S> /a 2 )(T c0 -T)/T is 
the mean-field perpendicular upper-critical field, and where 72 ~ (cos (^ + i)^ 2 -(/cbT c0 /-E/o) 1 / 2 | 
is an effective anisotropy parameter. The selective high-temperature expansion (55) indi- 
cates that this decoupling line is depressed quadratically by parallel field at temperatures 
outside of the 2D critical regime. Such behavior is consistent with anisotropic Ginzburg- 
Landau theory. 49 Continuity with the contour (60) in the coupled phase also indicates that 
the first-order line crosses the second-order melting line, T = T m . The phase boundary 
then continues up in field along the latter second-order line. Finally, Eq. (43) implies that 
the decoupling contour (59) continues into the coupled phase at temperatures T < T m 
down to the Josephson energy scale Tj. The non-zero line tension (42) for Josephson 
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vortices that is characteristic of the coupled phase results in a linear depression of this 
contour with parallel field [see Eq. (43) and ref. 50]. Last, it is possible that a vestige of 
the second-order melting transition found at high fields B± > B]_ persists down into the 
low-field region in the form of a crossover (see Fig. 1). This regime, however, is beyond 
the scope of the weak-coupling approach elaborated here. 

B. Latent Heat and Josephson Plasma Resonance. We shall now estimate the latent 
heat across the first-order decoupling line (59). The perpendicular contribution is equal to 
AE±/Af' = — J_ L A(cos(/)i ) i_|_i). The corresponding entropy jump is then equal to AS± = 
AE±/T. The entropy jump per vortex, al x AS±/Af'a 2 , at the first-order transition due to 
the Josephson coupling is therefore given by 

AS ± per vx = f3\\{B*J H D ){-A{cos^ l+1 ))k B . (63) 

Eq. (61) implies that it is of order — A(cos (fiij + i)kB- Last, it is perhaps useful to point 
out that a jump 

A{N[n z ])/N' = #lA(cos0 IiI+ i) (64) 

in the number of inter-layer fluxons per rung at the first-order decoupling transition is 
implied by the relation (14). 

Apart from playing an important role in the determination of the phase diagram 
discussed above, the local Josephson coupling between layers can also be directly probed 
through c-axis Josephson plasma resonance (JPR) experiments. Theory dictates that the 
plasma frequency is given by u p i = cuo(cos (pi^i) 1 ^ 2 , where ujq is its zero-field value. 10 ' 11 ' 46 
The local Josephson coupling is given by Eq. (55) in the vortex liquid phase, where 
£vx ~ a vx - This yields Koshelev's result 10 ' 11 

(upi/uo) 2 = (f <S>o/B ± a 2 )(J ± /2k B T) (65) 

for the corresponding JPR contours at constant frequency u p i <C ujq, where fo is weakly 
temperature dependent and of order unity (<7o). Consider next the coupled phase at per- 
pendicular fields in the weak-coupling regime a <C a vx <C A . The second term on the 
right-hand side of Eq. (43) for the local Josephson coupling is then dominant. The vortex 
contribution to the 2D correlation exponent (19) can be neglected in the low-temperature 
limit [see Appendix A, Eq. (A8)]. We therefore obtain the corresponding JPR contours 

(iv pl /u ) 2 = /o(5j/5±) T/T *° (66) 
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at constant frequency in the coupled phase at low temperatures, Tj <C T < Tm , where 
Bj is a constant perpendicular field scale of order B]_. Notice then that the perpendicular 
field along a given contour diverges exponentially with T* /T in the weak-coupling regime, 
uj p i <C ujq [see also Eq. (60)]. 

C. Bulk Limit. The thermodynamic limit, N — > oo, of an infinite number of layers 
must be taken in order to model extremely type-II layered superconductors in bulk. Since 
the decoupling contour (59) is determined primarily by the Josephson coupling in between 
adjacent layers, it should not change much in the case of an infinite number of layers. The 
first-order decoupling transition should therefore survive the bulk limit. 51 With respect to 
the question of the survival of the second order melting line in the bulk limit, it is useful to 
observe that the Coulomb gas ensemble (21) used here to describe the frustrated XY model 
also describes the unfrustrated model if one makes the replacement a vx — > a for the natural 
ultraviolet scale. Yet it is obvious that the universality class of the unfrustrated XY model 
should pass from 2D to 3D as N — > oo. Given that both the frustrated and unfrustrated 
layered XY models have a common Coulomb gas form (21), then the universality class of 
the former should pass from the 2D to the 3D universality class of the latter as N — > oo. 
In other words, the second-order melting line should survive the bulk limit as well. 

D. Bulk Pinning. Vortices both in real superconductors and in the XY model (1) 
can be pinned, respectively, by defects and by the underlying cubic lattice at low enough 
temperatures. To get an idea of how pinning can effect the Josephson coupling, we shall 
consider the optimal case where each c-axis flux line is fixed to some type of a correlated 
pin. The latter could represent a twin boundary, a columnar track, or a model substrate 
potential. If we now return to the analysis of the CG ensemble (21), this means (i) that 
the auto-correlation functions (16)-(18) for each layer in isolation are all identical, and 
(ii) that they resemble the auto-correlation function of the unfrustrated XY model on 
the square lattice due to the fact that vortex fluctuations are suppressed. In particular, 
they show quasi-long range order (£ vx = oo) at temperatures below the vortex/anti- vortex 
unbinding transition, ksTc 20 ^ = §</||, following Eq. (17). The phase autocorrelations 
within isolated layers are short range (£ vx < oo) at high temperatures T > Tc 2D \ on 
the other hand, following Eq. (18). The previous analysis holds in general, with the 
exception that the Coulomb gas ensemble (21) is not coarse grained. In particular, the 

24 



replacement a vx — * a must be made globally. This yields a field-independent result for 
the local Josephson coupling (13). No first-order phase transition is therefore expected. 
A second-order decoupling transition, however, remains at a critical temperature T c that 
should lie within the dimensional crossover window Tc 2D ^ < T < T x (see ref. 34). By 
the previous arguments, its universality class should pass from that of the 2D XY model 
without frustration to that of the 3D one in the thermodynamic limit of an infinite number 
of layers. A schematic phase diagram that summarizes the above results is shown in Fig. 
2. 

It must be emphasized that the above results are obtained under the assumption of 
optimum correlated pinning, which may not be realizable at all temperatures and fields. 
Indeed, an isolated 2D vortex lattice within the XY model is pinned only at low tem- 
peratures T < T^ D ^ inside of the quasi-ordered phase below the 2D melting temperature 
Tm D \ 29,30 The 2D depinning temperature, T^ D \ is null in the zero-field limit, and it 
increases linearly with increasing field. This depinning line, furthermore, merges with the 
2D melting line T = at a field of order $ /36a 2 . Since the focus of this paper is 

the mixed phase of type-II superconductors, which effectively have no substrate potential, 
we shall defer the analysis of intra-layer intrinsic pinning effects in the XY model (1) to 
future work. 



IV. Discussion 

Below, we shall compare the results obtained in the previous sections for the ther- 
modynamics of the mixed phase in extremely type-II layered superconductors with nu- 
merical simulations of the frustrated XY model, with relevant experimental work on 
high-temperature superconductors, and with the flux-line entanglement ideas developed 
by Nelson 52 and co-workers. 13 

A. Numerical Simulations. Let us first review the Monte Carlo simulation results that 
have been obtained recently for the uniformly frustrated XY model on an isolated square 
lattice. 29 Both these and simulations of the associated lattice Coulomb gas 30 find a vortex- 
lattice melting transition at a temperature (57) that is independent of field, in accord with 
theoretical expectations. 31 They also find, however, that pinning due to the square-lattice 
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grid at commensurate flux densities can only be neglected at large inter-vortex scales 



a vx > 6a. 



(67) 



(See the discussion at the end of section III.D.) Monte Carlo simulations of the layered 
XY model (1) have also been performed. 14 The most extensive simulations find evidence 
for a single first-order phase transition that coincides with the decoupling contour (59). 15 
The results reported in ref. 11, for example, yield a value for (cos^z+i) d near 1/2. These 
simulations find no evidence for a critical endpoint where the first-order transition termi- 
nates with decreasing temperature, however, contrary to the theoretical results obtained 
here (see Fig. 1). 

The cause for the above discrepancy may have two sources. First, the coarse-grained 
Coulomb-gas ensemble (21) used here to analyze the coupled phase at low temperatures 
may not be valid in regimes where the 2D vortex lattices are pinned by the XY model 
grid. This suggests that comparison of theory with the Monte Carlo simulations of the 
layered XY model (1) should be restricted to low perpendicular fields (67). Second, the 
phase diagram shown in Fig. 1 is valid only for anisotropy parameters inside of the range 



for XY layers of dimensions L by L. The lower bound guarantees the inequality Tj < 
Tm D \ while the upper bound insures that the system is not effectively decoupled at zero 
temperature. 34 Most, if not all of the Monte Carlo simulations performed to date have not 
satisfied conditions (67) and (68) simultaneously. This indicates that a comparison with 
the present theory is premature. 

B. High-T c Superconductivity. The phase diagram shown in Fig. 1 for the layered XY 
model is consistent with the nature of the flux-lattice melting observed in the mixed phase 
of clean high-temperature superconductors. 1-3 The latter melting transition begins at (or 
near) the zero field critical point, but it terminates in the middle of the T — H± plane as the 
temperature is lowered. The multi-critical point obtained here at To ~ Tm and Hq ~ B* ± 
is roughly consistent with this effect. It must be remarked, however, that the location of the 
critical endpoint observed experimentally in high-temperature superconductors is sensitive 
to the degree of disorder. 5 



5 < 7 < L/27ra 



(68) 
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The critical endpoint for the vortex-lattice melting line observed in clean high-temperature§ 
superconductors also coincides with a multi-critical point. In particular, the flux-line lat- 
tice frees itself from pins through thermal fluctuations across a vertical line in the T — H± 
plane that extends up from the critical endpoint. 5 ' 6 The combination of surface barriers 
plus the second-order vortex-lattice melting transition at T = T m that is obtained here 
theoretically for extremely type-II layered superconductors could account for this effect 
(see Fig. 1 and ref. 1). 

A so-called second peak line also extends horizontally to lower temperatures in the 
T — H± plane from the critical endpoint in clean high-temperature superconductors. 7 Bulk 
pinning is observed to increase dramatically at perpendicular fields that lie above this 
line. Recent observations of muon spin resonance (juSR) in the same systems have reached 
similar conclusions. 53 It is conceivable that this effect is related to the dimensional cross- 
over transition that is predicted here in the coupled phase between a flux-line lattice phase 
at low fields and a smectic solid phase at high fields (see Fig. 1). In particular, point 
pinning is relatively more effective in the latter phase composed of weakly coupled 2D 
vortex lattices. Note, however, that the cross-over predicted here for a finite number of 
layers in the absence of pins is possibly broad [see Eq. (60)], whereas the second peak 
transition is sharp. 5-7 This suggests either that a true decoupling phase transition occurs 
in bulk (see ref. 51), or that point pins are what drive the second-peak phenomenon into 
a sharp transition. 53 

C. Topological Defects and Flux-line Entanglement. Distinct solid phases of vortex 
matter have been proposed in the past. 4 In particular, Frey, Nelson and Fisher have pre- 
sented arguments for the existence of a "supersolid" phase in the mixed phase of layered 
superconductors that is characterized by highly entangled flux lines that preserve the 2D 
lattice structure within layers. 13 The "supersolid" phase is predicted to occurs whenever the 
temperature, Td, above which interstitial/vacancy pairs proliferate 12 lies below the melting 
temperature, T m , of the quasi-2D vortex lattice. The phase composed of Josephson-coupled 
2D vortex lattices that appears at low temperatures and high fields in the phase diagram 
for the layered XY model displayed in Fig. 1 evidently has physical characteristics that co- 
incide with the "supersolid" phase. 51 In particular, the flux arrangement is ordered within 
planes yet disordered along the field direction. The quasi-2D solid phase discussed here is 
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furthermore limited to temperatures inside of the window T D < T < T m at relatively high 
fields H± 3> Bj_. If the (Josephson) decoupling temperature scale T D is then identified 
with the scale Td, at which point vacancy/interstitial pairs begin to proliferate, 12 we re- 
cover precisely the "super solid" regime that is discussed in ref. 13. Last, we remark that 
an interstitial/vacancy pair in a single layer is confined by Josephson vortices that run 
between the neighboring layers. 12 This indicates that J± is the limiting energy scale for 
such defects. This in turn makes the identification of T D with T d plausible, since ksT D is 
also of order the Josephson scale J± in the high-field limit by Eqs. (43) and (59). 

V. Concluding Remarks 

A schematic phase diagram (Fig. 1) for the mixed phase of extremely type-II layered 
superconductors has been proposed on the basis of a weak-coupling duality analysis of the 
layered XY model with uniform frustration (1). The artificial pinning effects that arise 
from the 2D model grid in each layer were effectively removed by a coarse- graining (21) 
of the resulting dual Coulomb gas description (12). We notably predict that a smectic 
(super 13 ) solid type of vortex matter exists at high fields and low temperatures, and that 
this phase undergoes a second-order melting transition into a vortex liquid phase at higher 
temperatures. The solid phase exhibits a macroscopic Josephson effect, while the liquid 
phase does not. The weak-coupling analysis thus indicates that neither an intermediate 
"line-liquid" phase that exhibits a macroscopic Josephson effect between normal layers, 
nor a Friedel scenario characterized by decoupled superconducting layers, are thermody- 
namically possible in pure extremely type-II superconductors. It has also been argued that 
the second-order melting transition converts itself into a first-order one at lower fields of 
order many times the dimensional cross-over scale, B]_. This potentially accounts for the 
critical endpoint for the first-order melting of the flux-line lattice that is observed in clean 
high-temperature superconductors. 1-3 

Notably absent in this theory are effects due to magnetic screening 54 and to disorder 
caused by point pins. 55 Both of these effects are present to some degree in the mixed phase 
of all high-temperature superconductors. 5 ' 7 It remains to be seen how such departures 
from the pure XY model description (1) for the mixed phase of extremely type-II layered 
superconductors may change the conclusions reached above. 
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Appendix A: Villain Model 



The duality analysis for the conventional XY model presented in section II can be 
repeated in its entirety for the case the XY model in the Villain form: 



4 3) [?]= E n ^ 



WO)} 



^2 ^v n v\r -p(l 



exp 



2/?|| U + 2(3. 



(Al) 



In the weak-coupling limit, j3±_ — > 0, partial resummation of the type outlined at the 
beginning of section II. A yields the factorized form 



4 3) b] = ^cGb]-n^4 2) [o] 



for the partition function in terms of the corresponding 2D partition functions 



(A2) 
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2/9, 



1 X 

— ^n 2 (f,/) -i^2n(fj) ■ A{r) 



(A3) 



and the Coulomb gas ensemble (12). In contrast with the case of the conventional XY 
model, the fugacity is now equal to 



while the 2D auto-correlation functions (10) are set by the quotients 



(A4) 



Ci[qi] = Z^[ qi }/Z^[0] 



(2)r 



(A5) 



Repeating the derivation of the renormalized LD model (26) results in a renormalized 
Josephson penetration length (27) that displays an artificial temperature dependence ac- 
quired from the fugacity (A4). Also, Koshelev's formula (51) indicates that the local 
Josephson coupling acquires the same artificial temperature dependence in the decoupled 
phase. 

The above results were initially derived in ref. 19 by following an alternate route, 
where the spin-wave and the vortex degrees of freedom within layers are explicit. The 
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Villain model displays a complete factorization of such degrees of freedom. This permits 
the 2D auto-correlation functions (A5) to be computed analytically in the low-temperature 
ordered phase (see refs. 17 and 18). For example, the two-point auto-correlation function 
(15) for an isolated layer / has the asymptotic form 

e -*J>»^ (A6) 

in the ordered phase as a function of the separation r\i = \r\ — r^l, where A 1 is a suitably 
gauge-transformed vector potential, where r]2D = V + lvx is the 2D correlation exponent, 
and where ro is of order the natural ultra-violet length scale a vx set by the spacing between 
vortices. The latter is in accord with general scaling considerations. 37 The separate spin- 
wave and vortex components of rj2D are given by 

Vsw = (2n(3 l] )- 1 (A7) 

and by 

»?vx = ~ £ |rf <<JQz(0)<$Q,(f)>vx, (A8) 

r=£0 

respectively, where 5Qi(r) = Qi(f) — Qo(r) is the deviation of the intra-layer vorticity 
Qi(r) at the dual site r in layer / with respect to the flux-line lattice, Qo(r), that threads 
each and every layer at zero temperature. (Surface barriers are assumed.) The result 
(A8) is obtained through a cummulant expansion of the dipole pairs so indicated over a 
2D Coulomb gas (vx) with a uniform background charge density (Qi)/a 2 = b±/2na and a 
temperature of rj sw . Last, the prefactor in expression (A6) is of order 

<7o~(a/a vx )^. (A9) 

The upper bound 7q < e 27r ^n on the anisotropy parameter guarantees that go be of order 
unity for perpendicular fields B± > B]_. 

The methods used to compute the above two-point function (A6) in the ordered phase 
can be extended to compute the asymptotic form of any general n-point autocorrelation 
function (10). One thereby obtains the result 



|afe]| = <7o + ( °exp 



2 ^2 V2D- gj(ri) ln(|fi - f 2 \/r )qi(f 2 ) 

ri ^r 2 
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(A10) 



where n + (l) denotes the number of positive (or negative) probe charges applied to layer I. 
Notice that this expression is equal to the product of two factors, each corresponding to the 
spin- wave and to the vortex components of the 2D correlation exponent r]2D- The vortex 
factor in the generalized auto-correlation function (10) is obtained through the cummulant 
expansion of the associated Coulomb gas ensemble mentioned above [see Eq. (A8) and 
refs. 17 and 18]. 

With respect to the evaluation of the 2D auto-correlation function in the disordered 
phase at high temperatures, it is more useful to study the original dual form (A3) for the 

(2) 

corresponding partition function, (ri,r 2 )- The image of a string with its ends tied to 
the points r*i and r*2 at a low dual temperature (3\\ indicates the form 

4 2) (ri,f 2 ) oc e -<W0|i e - i /i 2 * f >- di? 

for the latter, where cry denotes the string tension (see refs. 35, 38 and 39). We thereby 
obtain the asymptotic result 

for the auto-correlation function of the 2D Villain model in the disordered phase, where 
£vx = /3||/cr|| is the correlation length. The string image can also be used to compute general 
n-point autocorrelation functions (10) in the limit where the probe charges correspond to 
well separated dipoles. The auto correlator in such case is simply given by the product of 
the two-point functions (All) represented by each dipole. 
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Appendix B: Fermion Analogy 



We shall now determine the thermodynamic character of the renormalized LD model 
[Eqs. (25)-(27)] in parallel magnetic field -B^+i = (^o /2tv d)bi t i + i in between layers / and 
I + 1. This model is defined by the energy functional 

r f N 1 7V_1 1 

£ld=J|| y dV|^-(V^ 2 -A - 2 ^cos[^ +1 (r)-^(r) + ^+^]| (£1) 

of the real phase variable Oi(r). It is known to be equivalent to coupled chains of spinless 
fermions at zero temperature, where each chain corresponds to a layer. 41 ' 42 (Previous analo- 
gies for layered superconductors in parallel field employed fermions that live in between 
consecutive layers. 56 ' 57 ) Specifically, the Hamiltonian for the fermion model is divided into 
intra-chain and inter-chain pieces, H = H\\ + H±, with respectively 
n r 

f dx v ¥ (*\id x * L - viidxVn) + U^[^ L ^ R (B2) 



and 



i=i 



N-l 



--U ± / dx 



i=i 



*l(x, l)^ R {x, I + l)* L (x, / + l)* R (x, I) + H.c. 



(S3) 



Above, ^ r(x, I) and *&l(x, I) represent the field operators for right (R) and left (L) moving 
fermions. The coordinate along the Josephson vortices, y, is related to the imaginary time 
variable r of the fermion analogy by y = v' F r. Here, the Fermi velocity v' F = v F sech2</> is 
renormalized by the intra-chain interaction t^, 58 with tanh2(/> = U\\/2txvy- Also, U± > 
is a repulsive backscattering interaction energy 58 in between chains. The Gibbs free-energy 
of the LD model (Bl) with respect to the normal state is then found to be related to the 
ground-state energy E-p of the fermion analogy by 



The identifications 



G cond /k B T = (L y /v' F )[E F (U ± ) - E F (0)}. 

b ljl+1 =2ir(N l+1 -Ni)/L x , 
V2D =2e 2 t 

A o 2 =ao 2 (\U±\/nv F )ri2D, 
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(B4) 

(55) 
(56) 
(B7) 



complete the equivalence between the models, where Ni/L x gives the fermion density in 
the Z th chain, where r]2D = ksT l2nJ\\ is the 2D correlation exponent, and where a is 
the natural ultraviolet length scale of the LD model (Bl). [Note that the above fermion 
analogy can directly be shown to be equivalent to the CG ensemble (21) in the screened 
phase, £ vx = oo, without recourse to the renormalized LD model (Bl) (see ref. 41).] 

Consider first the minimal case of iV = 2 layers. The above analogy then reduces to 
a Luther-Emery (LE) model for pseudo spin-1/2 fermions, 58 which exhibits pseudo spin- 
charge separation. The right-moving and left-moving spinless fermions corresponding to 

the pseudo-spin sector, ajjO) and b\\0), are governed by an ideal Hamiltonian along the 

t t t t 

Luther-Emery line t) 2 d = 1 of the form H a = Y^k[ v FH a k a k ~ K b k) + A a a ] k b k + A*6[a fc ]. 

They exhibit a spectrum ef(k) = ±(v' F 2 k 2 + A 2 ) 1 / 2 that is characterized by a pseudo spin 

gap 

A a = U ± /27ra . (£8) 

In the absence of parallel field, 6i )2 , the hole band (— ) is filled and the particle band (+) 
is empty. Since the pseudo-charge sector is independent of the coupling between chains, 
U±, the groundstate energy density relative to the decoupled state, U± = 0, is therefore 
equal to 

[E F (U ± ) - E F (0)]/L x = -A 2 j4nv F . (£9) 

(The logarithmically divergent contribution above that is equal to L" 1 J2k^<^( a t^k) has 
been removed by a normal-ordering type of procedure.) By Eq. (B4), the condensation 
energy per vertical rung at the corresponding LE temperature, /cbTle = 27rJ||, is then 
equal to 

-GJLCZLe)/^ = (aK/v' F f(J\\/2) (BIO) 

in zero parallel field. Yet in the limit of vanishingly small parallel field, the Gibbs free 
energy takes the form given by Eq. (36), where e\\ denotes the line tension for a single 
Josephson vortex. It is known that single-particle excitations within chains inherit the 
pseudo spin gap (B8) because of spin-charge separation. In addition, Eq. (B5) indicates 
that adding or removing a fermion from one of the two chains is equivalent to injecting 
a Josephson vortex in between the two superconducting layers. The equivalence formula 
(B4) thus yields the identity 

e ll (T LE ) = (\A rT \/v F )k B T LE (£11) 
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in between the line-tension and the pseudo spin gap. Comparison with Eqs. (B7) and (B8) 
yields the final expressions 

-^(Tle)/^' = (aa /A 2 ) 2 (J||/8) (512) 

and 

£||(Tle) =nJ ll a /A 2 (513) 

for the above free energies in terms of the original model (Bl) parameters. 

At arbitrary number of layers and at arbitrary temperatures, the fermion analogy can 
be treated in the mean-field approximation defined by the charge- density wave (CDW) 
order parameter Xl( x ) = (®\i( x i 1)^l( x , I)) and the associated gap equation 42 A; = U\\xi + 
U±(xi+i + Xi-i)- Standard self-consistent calculation yields a single-particle gap 

|Aj| = v F a^ 1 /smh[2nv F /(2U± - U\\)] = A (514) 

for all chains / that are not at a boundary (I = 1 or N), and for couplings that satisfy 
U± > U\\/2. Here, a\ is the ultraviolet cut-off for the mean-field theory, which is of order 
cko- At temperatures T > T* = A^J^jhs (or U\\ > 0), we have a null single-particle gap A; 
at small couplings U±_ < U\\/2 within mean field. Given the mean-field excitation spectrum 
s^(k) = ±(v F k 2 + | A;) 2 ) 1 / 2 for the fermionic excitations in chain /, the groundstate energy 
density relative to the decoupled state, U± = 0, is then equal to 

[Ef(U±) - E F (0)]/L x = -A 2 /4tiv f . (515) 

By Eq. (B4), the condensation energy per vertical rung is equal to 

-G ( l d /M' = (aA /v F ) 2 (k B T/4n). (516) 

for temperatures T > T*q in the vicinity of the decoupling transition. Eq. (B14) indicates 
that this quantity vanishes exponentially at the decoupling temperature T* that is set by 
the relationship U\\ = 2U±. 42 > 58 
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Figure Captions 



Fig. 1. Displayed is a schematic phase diagram for a finite number of weakly coupled XY 
layers with uniform frustration [Eq. (1)]. All artificial pinning effects due to the 
2D model grids are suppressed (see refs. 29 and 30). The decoupling (dashed/solid) 
line is a few to many times the dimensional crossover scale $ /7o 2a2 & t temperatures 
in the vicinity of T^n ^ [see Eqs. (60) and (61)]. Also, the Josephson temperature 
Tj = J±/k B is assumed to be smaller than the scale of the figure. Last, the solid 
phase should be defective on scales that are small in comparison to the Josephson 
penetration length Ao at elevated temperatures T > Tm . 

Fig. 2. Shown is a schematic phase diagram for the layered XY model (1) under the assump- 
tion that all of the vortex lines are fixed to some type of correlated pin. Vortex/anti- 
vortex pairs should begin to proliferate at elevated temperatures Tc 2D ^ < T < T c 
inside of the ordered phase on scales that are small in comparison to the Josephson 
penetration length Aq. 
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